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ABSTRACT 

Almost all of the > 600 known Kuiper belt objects (KBOs) have been discovered 
within 50 AU of the Sun. One possible explanation for the observed lack of KBOs beyond 
50 AU is that the distant Kuiper belt is dynamically very cold, and thus thin enough 
on the sky to have slipped between previous deep survey fields. We have completed a 
survey designed to search for a dynamically cold distant Kuiper belt near the invariable 
plane of the Solar system. In 2.3 deg^ we have discovered a total of 33 KBOs and 1 
Centaur, but no objects in circular orbits beyond 50 AU. We find that we can exclude 
at 95% CL the existence of a distant disk inclined by i < 1° to the invariable plane and 
containing more than 1.2 times as many D > 185 km KBOs between 50 and 60 AU as 
the observed inner Kuiper belt, if the distant disk is thinner than a = 1?75. 

Subject headings: Kuiper Belt — minor planets — solar system: formation 

1. Introduction 

More than 500 Kuiper belt Objects (KBOs) have been discovered since the first, 1992 QBl, was 
detected almost a decade ago. Nearly all of these are currently within 50 AU of the sun, and only 
one has been detected beyond 54 AU, despite surveys sensitive enough to detect 160 km-diameter 
KBOs to 60 AU. Such surveys have placed strong upper limits on the density of KBOs in an outer 
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Kuiper belt (Allen, Bernstein & Malhotra 2001; Trujillo et al. 2001; Jewitt, Luu, & Trujillo 1998; 

Gladman et al. 2001). There are several possible explanations for the observed absence of distant 
KBOs, each with distinctly different predictions for the Kuiper belt beyond 50 AU, and distinctly 
different implications for the history of the solar system. 

In the simplest view of the Kuiper Belt, the region beyond 50 AU, where the gravitational 
perturbations of the giant planets become negligible, has remained dynamically cold, with KBOs 
moving on their primordial circular and coplanar orbits. It is inferred that the region inside 48 AU 
has been heavily depopulated, by a factor of 10-100 in mass, over the age of the solar system as a 
result of planetary perturbations and mutual collisions (Stern 1996; Stern Sz Colwell 1997; Levison 
& Duncan 1993; Duncan, Levison, & Budd 1995). In the absence of such perturbations, the outer 
disk is expected to be a factor of 5-100 times as dense as the presently observed inner Kuiper belt. 
However, as yet no discovered KBO has been determined to be in a circular orbit beyond 48 AU. 
Shallow surveys, with typical magnitude limits of i? < 23.5, have searched for KBOs in more than 
100 deg^ of sky, but objects at > 50 AU must be > 280 km in diameter to be detected in such 
surveys. (In quoting sizes of KBOs in this paper, we assume a 4% albedo.) It is desirable to search 
for smaller objects since the growth timcscalcs for KBOs in the outer Kuiper belt, though highly 
uncertain, are likely longer than in the inner Kuiper belt. Only a few square degrees, however, have 
been surveyed to magnitudes R > 24.5. 

Hahn (2000) has pointed out that if the outer disk were dynamically very cold, thus restricted 
to a thin line on the sky, it could have easily slipped between the previous sparse deep survey fields. 
The most likely location for this distant cold disk would be the invariable plane, the plane normal to 
the total angular momentum of the solar system, which is inclined 1?6 to the ecliptic (Allen 1976). 
Most deep KBO surveys, including our previously published results (Allen, Bernstein &; Malhotra 
2001) [ABM], have been conducted close to the ecliptic plane, possibly missing a thin disk near the 
invariable plane. We have therefore supplemented the ABM survey with four new fields that are 
strategically chosen on the sky to maximize the likelihood of detecting a thin trans-Neptunian disk 
near the invariable plane. We present here new limits on a putative cold disk beyond 50 AU that 
result from our observations in these new survey fields combined with previously published data. 

2. Distant Disk Survey Observations 

The new observations were taken with the KPNO Mayall 4-m telescope and Mosaic camera, 
during two consecutive nights in March 2001. The observation and rcdTiction methodologies are 
identical to those of ABM. In brief: data were taken using a broad-band VR filter similar to that 
described in Jewitt, Luu, & Chen (1996). Each field is observed in a series of 7 to 9 exposures each 
night, each exposure 480 s long. The images were registered and then combined with a "digital 
tracking" technique to search for faint moving objects: after registration, the images from each night 
are summed (with sigma-clipping rejection) into a deep template. This template is then subtracted 
from each individual image in the opposing night, after matching the point spread functions. This 
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removes all stationary objects from each image, leaving only moving objects and cosmic rays. These 
blank images are searched individually for bright KBOs, then combined at potential KBO motion 
vectors to search for faint slow- moving objects. The images were combined tens of thousands of 
times to search all likely orbits between 30 and 100 AU. 

As the observations were taken in fields very near opposition, the distance can be estimated 
to within 10% from the apparent velocity. Although we cannot determine orbits with the short 
(two-night) observational arc obtained here, a distance alone is usually sufficient to exclude the 
possibility of a distant circular orbit. 

We discovered 10 new KBOs in this 2001 survey data, over an effective area of 1.0 deg^. 
The effective area is simply the subtended survey area multiplied by the detection efficiency — 
the latter factor being probability of detecting a KBO at a particular magnitude. The effective 
area is determined by inserting hundreds of artificial moving objects into the data, then repeating 
the search process. In our 2001 survey, the detection efficiency drops by 50% at a magnitude of 
R ~ 24.9 due to noise, but has no dependence on velocity or distance as all KBOs within 200 AU 
will move appreciably in the two-night baseline. This limiting magnitude varies from field to field 
due to seeing and integration-time variations, while the peak detection efficiency varies primarily 
with stellar density and field overlap. Data for each field are listed in Table 1 and the details for 
each discovered object are listed in Table 2. 

Also included in the Tables are details of the fields and detected objects from the 1998 and 
1999 ABM data, which have an effective area of 1.3 square degrees with mean limiting magnitude 
of i? = 25.4. Altogether, our surveys between 1998 and 2001 cover 2.3 square degrees and we 
discovered a total of 33 KBOs and 1 Centaur. The discovered objects were found at distances 
between 20 and 53 AU, but none of the objects beyond 48 AU are on circular orbits. In our 2001 
survey, no KBOs were found beyond 46 AU. We have clearly detected no possible distant cold disk 
members, even though we could have detected KBOs larger than D=185 km to distances of 55 AU 
or greater over nearly the entire survey area with our stated magnitude limits. 

Figure 1 shows in invariable-plane coordinates the extent of the ABM fields (A-G), the new 
fields (H-L), and other authors' surveys. The heavy line is the ecliptic plane (the plane containing 
the orbit of the Earth), about which the previous surveys are clustered. Observing in the ecliptic 
rather than the invariable plane (the total angular momentum plane of the Solar System) has been 
the norm, as the difference between the two is not important when considering a thick disk such as 
the Kuiper belt inside 50 AU. It is only when searching for a disk with vertical dispersion a < 2° 
that this diflFerence becomes significant. 

Deep KBO surveys conducted by other researchers (Chiang &; Brown 1999; Gladman et al. 
2001, 1998) generally cover too little area (on the order of a Keck field ~ 0.01 deg^) to provide 
additional constraints on a thin disk. In the analysis described in the next section, we have included 
the KBO survey by Gladman et al. (2001), conducted with the Canada- France-Hawaii Telescope. 
This survey discovered 1999 DG8, for which a one-night arc suggests a distance of « 60 AU. This 
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is the only known object which could have a circular orbit beyond 50 AU, but the short arc means 
that the eccentricity of the orbit is completely indeterminate — it could easily be a scattered disk 
member. For the following analysis we must calculate the probability of cold disk models being 
consistent with the observations. We sum the probabilities of the model yielding 1 or objects in 
the Gladman et al. (2001) CFHT field, since the nature of 1999 DG8 is indeterminate. This, along 
with the high latitude of the field [1.5 to 2 deg to the invariable plane], substantially weaken its 
ability to constrain a cold disk in the invariable plane. 



3. Limits on a Distant Disk 

We did not detect any KBOs beyond 50 AU in the primordial, circular orbits expected of 
objects in a distant cold disk. Wc wish to constrain the total number of KBOs which could be 
resident in such a disk. We will consider only KBOs with diameter D > 185 km, since such objects 
are detectable at > 50% efficiency in all our fields at a distance of 55 AU, assuming 4% albedo. We 
express the total population N^oid of D > 185 km objects between 50 and 60 AU as /-/VckbOj where 
-^VcKBO is the number of classical KBOs (CKBOs) larger than D = 185 km. Trujillo et al. (2001) 
determine A'^cKBO ~ 3500. The cold disk is assumed to have some inclination i and ascending node 
0, relative to the invariable plane. The surface density of cold disk objects on the sky is assumed 
to have a Gaussian distribution in the vertical direction with dispersion a, and to be uniform with 
disk longitude. An exponential vertical distribution yields similar results for most disk widths. 

Wc ask then what is the probability P{0\f,i,Q, a) of a given disk model yielding the obser- 
vations O. For simplicity we assume that all the A'coid objects in the 50-60 AU range have the 
apparent magnitude m = 24.8 that a D = 185 km object would have at a distance of 55 AU. This 
is conservative in the sense that many of the objects would presumably be larger and brighter, but 
we do not want to assume any particular size distribution for the cold-disk members. Note that 
for / = 1, we may imagine taking the CKBO population, collapsing (cooling) it down to a thin 
disk, then moving the disk out from its present 38-48 AU annulus to a 50-60 AU annulus. If the 
"classical" belt is truly a depleted and excited inner portion of the full Kuiper Belt, then this is a 
lower bound to the expected distant disk. 

For chosen disk parameters {/, i, f2, a}, we convolve the model sky surface density with the 
survey geometry and detection efficiencies specified in Table 1, yielding an expected number of 
detections. The expected number of detections for each field, Ni, for each cold disk orientation is 
then 

(1) 



AT, =3550^^^^^^^°^^^^^) 



2 27r 



where Ei is the detection efficiency in each field, A(RA) is the field width and bn and bi2 are the 
top and bottom field latitudes in the cold disk plane. Since the actual number of detections is zero 
(save possibly for 1999 DG8 in the CFHT field), P{0) is easily calculated from Poisson statistics. 
A given disk configuration is excluded at 95% confidence if > 3 detections are expected. 
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As an example, 74% of the nodes of a very thin disk {a = 0?5) inchned at i = 1?0 to the 

invariable plane, are ruled out at the > 95% confidence level for / = 1. Most of the remaining 
un-excluded disks travel through our field H, due to its slightly brighter limiting magnitude. The 
i = 1°, (T = 0?5 disks which arc not ruled out by our deep fields are plotted in Figure 1. 

To simplify the interpretation, we marginalize over the ascending node by assuming a uniform 
prior distribution in the interval < < 27r: 

1 f"^^ 



Figure 2 plots the probability P{0\f = 1, i, a). As we allow the cold disk model to be more inclined 
to the invariable plane, our limits become weaker because these disks more easily dodge the surveyed 
fields. Our limits also become weaker as the model disks become thicker, at inclinations < 1°, since 
the sky density of the KBOs drops in our fields. This situation reverses at higher inclinations, 
where the thicker disks (cr > 0.5°) are more likely to be detected, as they are unable to slip between 
our fields and the inclination variation matters less. 

If we assume that the distant cold disk must be located in the invariable plane, then we can 
place further limits on how many objects larger than 185 km could exist between 50-60 AU without 
being detected in our survey. Figure 3 plots the value of / at which P{0\f,i = 0°,a) = 0.05, as a 
function of disk thickness a. A thin disk (a = 0.22°) must have less than half as many D > 185 km 
KBOs as the classical Kuiper belt (/ < 0.5 at 95% CL). An i = 0° distant disk with the same 
number of objects as the inner belt must have a > 1.85° to have a > 5% chance of escaping 
detection in the surveys to date. The apparent width of the classical Kuiper belt is ^ 5° (Trujillo 
et al. 2001). If the outer disk has the same width instead of being dynamically colder, an increase of 
factor / > 2.36 over the classical belt population is ruled out at the 95% CL. The case of equivalent 
inner- and outer-disk sky distributions in considered in more detail in ABM. 

By also marginalizing over the inclination, we can place similar limits on / for cold disks within 
1° of the invariable plane. Assuming that the model disk pole locations are uniformly distributed 
in polar area, then 



We can then find / at which P{0\f, i < 1, a) = 0.05, plotted in Figure 3 as the solid line. We can 
place tighter limits on a disk in the invariable plane, as above, than on any disk in any inclination 
up to 1° because more of our fields were located in the invariable plane. However, we find for 
a < 1.75° we can exclude at 95% CL or better the existence of a distant disk at any inclination 
i < 1° to the invariable plane that contains at least 1.2 times as many D > 185 km KBOs as the 
Classical Kuiper belt. 




(2) 



P{0\f,i<l°,a) 



£ P{0\f,a,i)idi 



(3) 
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4. Conclusions 

We have provided strong limits on the existence of a distant cold disk. The simple expectation 
of a dense thin disk in the invariable plane appears to be ruled out. However, it is possible that 
a thin disk composed of objects smaller than D = 185 km, undetectable in much of our survey, 
is present. This would rcqTiirc an extreme change in the size distribution of KBOs jiist beyond 
~50 AU, as objects larger than D ^ 1000 km have been found at distances near 50 AU. It is not 
clear what would cause such a large and sudden change in the size distribution. 

Other explanations for the observed lack of distant KBOs involve dropping the density of 
objects beyond 50 AU. One possibility is that the outer Kuiper belt was dynamically excited by a 
stellar encounter early in the history of the solar system (Ida et al. 2000). An increase in the mean 
eccentricity and inclination of distant KBOs, thus lowering their apparent sky density and possibly 
halting accretion of large KBOs, would make a distant disk much harder to detect. Dynamical 
excitation by a stellar encounter results in a distinctive orbital distribution of KBOs, and could 
yield limits on the birth cluster environment of the solar system (Adams k. Laughlin 2001). 

It is also possible to explain the lack of KBOs beyond 50 AU if the primordial planetesimal disk 
ended at this distance. If Neptune and Uranus have migrated significantly over the age of the solar 
system (Malhotra 1993, 1995; Thommes, Duncan, & Levison 1999), the primordial solar nebula 
surface density beyond ~ 30 AU could have been very small initially. This could suggest that the 
larger KBOs actually formed interior to 30-40 AU and were displaced to their present orbits by 
a large scale rearrangement of orbits in the outer Solar system. Circumstellar disks around other 
stars have been observed to have diameters ranging between 50-1000 AU (McCaughrean k, O'Dell 
1996; Brown et al. 2000). A truncation of the solar nebula near 50 AU would fit easily within this 
range, with important implications for planetary formation models of our Solar System. 

In order to distinguish between these theories and further constrain the distant KB population, 
deeper and wider survey observations of the Kuiper belt (to limiting magnitudes of i? > 25) will 
be necessary. Based on the limits we have found here, a survey covering approximately 20 square 
degrees near the invariable plane should be able to detect a distant disk dynamically excited to 
a thickness of 10° with even half the number of objects as the CKB. This would constrain the 
density and distribution of the Kuiper belt at distances beyond 50 AU, including the Scattered 
Disk Objects, enabling a better comparison between our solar system and extrasolar systems. 

Wc thank L. Strolger and M. Holman for attempting to recover some of these objects, and 
the staff at Kitt Peak National Observatory for their excellent support. This work is supported by 
NASA Planetary Astronomy grant #NAG5-7860; GB is further supported by grant #AST-9624592 
from the National Science Foundation; RM is further supported by NASA grants #NAG5-10346 
and #NAG5-11661. 
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Invariable plane longitudE (degrees] 



Fig. 1. — Fields from this survey, and others deep enough to provide hmits on the location of distant 
thin disk aX R = 24.8, are shown here in invariable plane coordinates. The heavy line shows the 
location of the ecliptic plane. The lighter curves show the locations where a distant disk, inclined to 
the invariable plane by 1°, might escape detection in surveys to date. The distant disk is assumed 
to have a Gaussian thickness of 0?5, contain as many D > 185 km objects as the "classical" Kuiper 
belt, and be excluded if there is < 5% chance of being consistent with our observations. 
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Fig. 2. — The probability of detecting a distant disk with varying thickness and incHnation (aver- 
aged over ascending nodes) to the invariable plane. This assumes that the number of objects larger 
than D = 185 km is similar to the number oi D > 185 km KBOs in the classical Kuiper belt. 
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Fig. 3. — Assuming the distant cold disk is in the invariable plane (i = 0°), we calculate the upper 
limit to the ratio / of D > 185 km objects relative to the classical Kuiper Belt, shown by the 
dashed line. By marginalizing over the inclination (averaged over inclination pole angular area) 
we can also calculate the upper limit to the ratio / for all cold disks with i < 1° to the invariable 
plane, shown by the solid line. Disk populations above each line are inconsistent at > 95% CL with 
our non-detection. 
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^Observations marked with asterisks were searched for KBOs; unmarked observations were used only for (p)recovery. 

and Eq are fit to the detection efficiency E for detection of implanted KBOs using E = -^crfc[(m — m^Q)/y/2w\. Eg is 
the peak detection efficiency for each field, m^o is the point where the efficiency is 50% of Eq, and ti; is a width parameter. 

'^Detection efficiencies for fields E and F are based on averages from nights 1 &: 2 and nights 3 & 4, as these fields search the 
same section of sky. 



Table 2. Objects Discovered 
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F 


9 d 


24.1 


32.4 zb 5.0 


r\ nn 1 n nn 

0.20 ± 0.22 


29.8 ± 6.4 


39.1 ± 1.4 


129 


8.2 


1999 JF132 


F 


9 d 


24.0 






1.6 ±0.4 


43.1 ±2.1 


164 


7.7 


1999 JH132 


F 


9 d 


25.5 






0.6 ±0.2 


41.1 ±2.4 


75 


9.4 


1999 JJ132 


F 


9 d 


25.5 






3.2 ±1.2 


50.1 ±2.8 


111 


8.5 


1999 JK132 


F 


9 d 


24.8 






16.0 ±5.9 


39.0 ±2.5 


93 


8.9 


1999 KT16 


F 


1 d 


25.1 






8.5 ±3.7 


46.3 ±3.1 


114 


8.4 


1999 KK17 


F 


9 d 


25.0 






8.7 ± 16. 


49.8 ± 7.5 


139 


8.0 


1999 KL17 


G 


90 d 


25.2 


46.2 ±0.2 


0.03 ± 0.17 


2.8 ±0.0 


47.6 ± 0.0 


115 


8.4 


1999 KR18 


G 


89 d 


24.9 


43.3 ±0.3 


0.21 ± 0.03 


0.6 ±0.0 


52.6 ±0.1 


162 


7.7 


G3-01^ 


G 


9 d 


25.6 


39.9 ±2.7 


0.16 ± 0.05 


1.6 ±0.1 


33.4 ± 1.4 


47 


10.4 



Table 2 — Continued 



MPC 


Field 


Arc 


ruR^ 




i 




Heliocentric 


Diameter^ 




Designation 




Length 




(AU) 


n 




Dist. (AU) 


(km) 




2001 FB185 


H 


1 d 


24.6 




4.5 ± 1 


5 


39.8 ±3.1 


106 


8.7 


2001 FC185 


J 


1 d 


23.8 




••• 21.2 ±8 


9 


39.0 ± 3.4 


147 


7.9 


2001 FE193 


K 


1 d 


23.2 




3.1 ±1 


1 


42.8 ±2.5 


234 


6.9 


2001 FG193 


K 


1 d 


24.2 




3.4 ±1 


6 


41.8 ±2.5 


141 


8.0 


2001 FD193 


K 


1 d 


24.2 




■•• 12.3 ±4 


6 


44.7 ±2.7 


161 


7.7 


2001 FC193 


K 


1 d 


24.8 




2.8 ±1 


5 


46.4 ±2.6 


132 


8.2 


2001 FJ193 


L 


1 d 


24.5 




1.1 ±0 


6 


36.4 ±2.4 


93 


8.9 


2001 FF193 


L 


1 d 


23.8 




2.0 ±1 


6 


44.7 ±2.6 


194 


7.3 


2001 FH193 


L 


1 d 


24.8 




8.9 ±3 


5 


40.9 ±2.6 


102 


8.7 


2001 FL193 


L 


1 d 


24.9 




1.0 ±0 


6 


40.9 ±2.5 


97 


8.8 



*A11 1999 and 2001 objects were detected in the VR filter. The listed magnitude may be transformed 
to standard RviaR = VR-0A6{V - R-0.5). As0.3<V-R< 0.7 for KBOS (Tegler & Romanishin 
2000), the R mag varies by 0.2 mag or less with color. 

'^No data are given when the arc is too short to provide meaningful constraint. 

^Diameters assume albedo of 0.04. 
'^Absolute magnitude in R 

^'Objects E2-01 and G3-01 were not reported to the MPC due to insufficient S/N on the recovery 
observations. 



